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We study the interaction of a Bose-Einstein condensate, which is confined in an 
optical lattice, with a largely detuned light field propagating through the condensate. 
If the condensate is in its ground state it acts as a periodic dielectric and gives rise to 
photonic band gaps at optical frequencies. The band structure of the combined system 
of condensed lattice-atoms and photons is studied by using the concept of polaritons. 
If elementary excitations of the condensate are present, they will produce defect states 
inside the photonic band gaps. The frequency of localized defect states is calculated 
using the Koster-Slater model. 



I. INTRODUCTION 

The achievement of Bose-Einstein condensation in magnetic traps [[!) has induced a great interest in 
the properties of quantum atomic gases and their manipulation by atom optic techniques. Although 
the latter are usually used for laser cooling in the formation process of a Bose-Einstein condensate 
(BEC), confinement in an optical dipole-trap has been demonstrated only recently 0. The all-optical 
confinement of a BEC provides a great potential for the manipulation and application of BEC. In 
particular, it opens the new opportunity to study atomic BECs in optical lattices. In recent years 
experimentalists have made great efforts to create a BEC in optical lattices. Although there are 
presently still some technical problems to achieve this goal the rich physics of uncondensed ultracold 
atoms in optical lattices || and quasi-crystals Q has attracted great interest for both experimentalists 
and theorists. Recently, several theoretical papers dealing with condensates in optical potentials have 
been published |||| 

In this paper we focus on a new aspect of this subject: light propagation through a coherent 
condensate confined in an optical lattice. Since the ground state of the condensate in a lattice potential 
is periodic, it will act as a periodic dielectric for laser light propagating through it. Thus it will give 
rise to photonic band gaps at optical frequencies. 

The phenomenon of photonic band gaps is a natural consequence of the periodicity of the condensate. 
In fact, it also should occur for uncondensed ultracold atoms in optical lattices. However, in the case 
of a condensed atomic lattice what is interesting is that, because of the macroscopic occupation of the 
ground state, a proper description of photonic band gaps is given in terms of polaritons (an entangled 
coherent system composed of superpositions of photons and excited atoms). Furthermore, elementary 
excitations may be present in the lattice BEC. In general these excitations are no longer periodic and 
will cause distortions of the perfect periodic structure of the condensed atomic lattice. An excitation- 
induced defect in the atomic lattice in turn causes the occurence of defect states inside photonic band 
gaps of light propagating through the BEC. In this sense, elementary excitations have a close analogy 
to lattice defects in solid state physics which cause defect states in the electronic band structure. 

The paper is organized as follows. In Sec. II we will derive the equations of motion. In Sec. Ill we 
consider the case where only the lattice laser beams are present and seek for a periodic solution to 
the coupled equations of motion describing the ground-state BEC and the lattice laser beams. This 
solution shows that the optical potential experiences no back- reaction from the condensate if the latter 
has settled down into its ground state. In this sense, the lattice laser beams just effectively act as 
a constant periodic potential for the BEC. In Sec. IV we consider the propagation of a weak probe 
laser beam through the ground-state BEC and derive the form of the lowest photonic band gap for 
this beam by using polariton modes. To examine the behaviour of a probe laser beam propagating 
through a weakly non-periodic BEC a theory of defect states for photonic band gaps is developed in 
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Sec. V which is applied in Sec. VI to the Koster-Slater model for a localized elementary excitation. 
Sec. VII concludes the paper. 



II. EQUATIONS OF MOTION 

The system under consideration consists of interacting two-level atoms which are coupled to the 
electromagnetic field. This coupling is described by using the electric-dipole and rotating-wave ap- 
proximation so that the corresponding second quantized Hamiltonian is given by 



H — Ha + i?NL + #E.M. + -Hint , (1) 



where 



H A := / d 3 x ]T *!{|^ + V(x) + Ei}% (2) 

i=e,g 

describes the atomic center-of-mass motion. V(x) denotes an external potential. M represents 
the atomic mass, and i — e,<? are the internal energy levels for ground-state and excited 
atoms, respectively. The corresponding field operators ^f g and f e fulfill the commutation relations 
[ty i(xy , *5> j (y)] — 5ijS(x — y). J/nl is the nonlinear part of the atomic Hamiltonian which describes 
two-body collisions. For a dilute Bose gas it can be approximated by 

Ha :=IJ d 3 x J2 . (3) 

where :— Awh 2 ai l J^ / M are coupling constants and denote the scattering lengths for scattering 
between atoms in the internal state i and j. 

For the description of the electromagnetic field we use the representation in terms of positive end 
negative frequency parts of the vector potential, A(x) = A^ + \x) + A^'(x). This representation will 
turn out to be convenient for the adiabatic elimination of excited atoms. The Hamiltonian for the 
free electromagnetic field then takes the simple form 

r 3 

H EM . = 2s / d 3 x £ A{-\uj 2 ) ab A { b +) . (4) 

J a, 6=1 

The positive and negative frequency parts are related by (A^Y — A^ and fulfill the commutation 

relation [A^ (x), A b (y)] — (h/2eo)oj~ 1 S a r b (x — y), where S^ b (x — y) is the transverse delta function. 
The frequency operator Co is a pseudo-differential operator whose action is defined in momentum space 

by 

[u>A](x) = (2vr)- 3 / 2 / d 3 ke lts c\k\A(k) . (5) 



The physical interpretation of the frequency operator is simple. It just multiplies a photon mode 
with its frequency u)(k) — c\k\. A more compact representation of cD in position space is given by 

Cj = cy/—A = c| — tV|, where A denotes the Laplace operator. 

Using the positive and negative frequency part of the vector potential the electric dipole coupling 
between the atoms and the electromagnetic field in rotating-wave approximation can be written as 

H int = i [d 3 x{*l* e (d* ■ CoA^) - *+* fl (d- ojA^)} . (6) 

The Heisenberg equations of motion derived from the Hamiltonian (Q) are given by 
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lk ^ e =ljKI + V + Ee+ T, 9ej^j \ *e - m 9 {d- U>&») (7) 



to*B = \ Jjj + V + E g + E to*]** f *5 + ' ( 8 ) 

[ J=e,g J 



6=1 



III. LATTICE LASER BEAMS AND BEC IN GROUND STATE: DECOUPLING OF THE 

FIELDS 

To analyse the interaction between the atoms and the lattice laser beams in absence of an exter- 
nal potential (V(x) = 0) we restrict us to the particular case where the atomic field is composed of 
condensed atoms, i.e., a Bose-Einstein condensate. This allows us to make further substantial sim- 
plifications. As is well-known a condensate can be described by assuming that all atoms are in the 
same quantum state ip g . This amounts in replacing the field operator ^f g in Eq. (^) by the c- number 
field ip g which then fulfills a nonlinear Schrodinger equation. In addition, we assume that the photon 
fluctuations of the lattice laser beams are small and therefore not important for our case. This allows 
us to replace the operator by a corresponding classical field vector A^ of the lattice laser beams. 

We consider the regime of coherent interaction where the electromagnetic field is detuned far away 
from the atomic resonance frequency uj rcB '■— [E e — E g )/%. Specifically, we assume that the detuning 
Al '■= lul—lUj-cs (with Ui := c|fcj,|) of the lattice laser beams is negative (red detuning) and its absolute 
value much larger than any other characeristic frequency of our system so that we can adiabatically 
eliminate the excited atoms [^-10|. This amounts in replacing the field operator for excited atoms by 



—i 



* e « ■ (10) 



Inserting Eq. (10) into Eqs. (||) and (||) one easily finds (to first order in 1/Al) 

f 
2M 



= {^ + V + E n + 9 9g ^ g + jr^id- CjA\ >)(«?• V ff (11) 



i{A^\x)) a = u>(A { L + \x)) a + I d 3 y\4> g (y)\ 2 (d- CuA ( L + \y))^2dtS^ b (S-y) (12) 
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Eqs.((T|) and (|T^) describe the coherent coupling of a ground-state atomic field to the lattice laser 
beams. The physics implicit in these equations is straightforward. The laser beams induce an optical 
potential for ground state atoms which is proportional to the light intensity (oc (ojAl) 2 ). The atoms 
in turn act on the electromagnetic field like a dielectric, where the index of refraction is determined 
by the density |V>g| 2 of ground-state atoms. 

We are interested in how the condensate affects the lattice laser beams and the corresponding back- 
reaction in the optical potential. For simplicity, we take the laser beam to be parallel to the x-axis. 
Further, we assume that the BEC has settled to its ground-state which, because of the periodicity 
of the optical potential provided by the lattice laser beams, is periodic. It is then convenient to 
decompose the fields into a discrete Fourier series 



^s( x ) = exp[ilkL ■ x] 

lez 

fi( £ >(ic) := ^d- CjA [ +\x) = ni i L) exp[ii*fc • x] 



lez 



We remark that Ul does not denote the wavevector of the laser beams. It is defined by its relation to 
the spatial period xl of the optical lattice by = e x 2ir /xl- This period Xi differs in general slightly 
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from the wavelength of the laser beams outside the atomic medium jy]. Transforming Eqs. (O) and 



(12) to momentum space we arrive at the following one-dimensional set of equations, 



^ = If ^ + tofc E + t^4a7 e (13) 

^ ' m,n£Z ^ ' m,«eZ 

.TiftW = ^lilflW + "$f\ E ^'C^i , (14) 

m,n£Z 

where we have defined the transversal dipole moment of the atoms d± := d — kr,(d ■ ki,)/k L . We 
can now exploit the fact that the optical frequency u>l is typically (very) much larger than any other 
frequency scale involved in the system. This allows us to perform a rotating-wave approximation by 
inserting Cl^ :— exp{ia;Li|Z|}f2j L ' ) into Eqs. ( |l3| ) and (jl^) and neglecting all terms which rotate at 
multiples of the frequency u>l- This procedure results in the simplified equations 




(2^) 3 A L 



^ + 2m^*^-l (15) 



mGZ 



i«if L) - huj L \i\ |i + —A^g w >} n w + _A_ WL | Z | n w |^ ; _ 2 ^ . (16) 



These equations have some properties which allow to decouple the system for some physically interest- 
ing cases. The most important one is that the atoms do only couple counterpropagating modes, i.e., 
n\ L) andfi^. This is a direct consequence of energy conservation since a transition to any other mode 
would require an amount of energy in the order of Tvusl which cannot be provided by the interaction 
with the condensate. In addition, it is not difficult to see that both the mean density p := ^ n \i[) n \ 2 

and the mean light intensity per mode, /; := | 2 + |fi[v| 2 , are conserved quantities, reflecting the 
conservation of the total number of atoms and the number of photons with energy respectively. 
Therefore, the first sums on the right-hand-side of Eqs. ( |l5| ) and ( pi| ) just produce a constant shift of 
the energy levels. 

We now consider a solution of the system of equations (|l5| ) and (|l6|) which corresponds to a standing- 
wave lattice laser beam interacting with a BEC in its ground state in the coherent regime. For a 
standing- wave lattice, we can make the ansatz f2, = ■ In addition, in the ground state of 
the BEC the time dependence of all coefficients i[>i 1S given by if>i(t) = exp[—ipt /h]ipi(0) , where fi 
denotes the chemical potential and ipi (0) can be chosen to be real. It is then not difficult to show that 
fij = holds for all times (by differentiating both sides and comparing the results). In addition, 

the expression Q{ L 'Q,_^*, which describes the optical potential in Eq. (|l5|), is time-independent, too. 

An immediate consequence of this fact is that in the coherent regime the optical potential created 
by the lattice laser beams for the BEC in Eq. (|l^) is not altered by the condensate itself. Thus, the 
problem is decoupled and the condensate behaves as if it were moving in a given external periodic 
potential. In this sense the influence of the lattice laser beams on the BEC can be replaced by an 
external potential V(x) with periodicity ir/kL- 



IV. POLARITON BAND THEORY FOR LIGHT INTERACTING WITH A CONDENSATE 



In the preceding section we have shown how the combined system of a BEC and lattice laser beams 
behaves in its ground state. We are now interested in a different situation where a running probe 
laser beam propagates through the "lattice condensate" (condensate plus lattice laser beams). The 
behaviour of the probe laser can intuitively understood by considering the BEC as a kind of dielectric 
for the probe laser beam. Since the BEC is periodic the probe laser beam will effectively propagate 
through a periodic dielectric. We thus expect it to show the phenomenon of photonic band gaps |l2| . 
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To describe the interaction of the probe laser beam with the condensate we will assume that the 
ground-state BEC changes little so that ijj g enters as a given external field into the equations of 
motion for excited atoms ^ and for the probe laser beam (||). Solutions of these coupled equation 
describe polariton modes, i.e., entangled superpositions of matter and light fields |D|[l4j . Thus, it is 
really "polaritonic" band gaps rather than photonic band gaps that we are studying. However, for 
sufficiently large detuning of the lattice laser beams the entanglement is very small so that the result 
indeed can be considered as photonic band gaps. 

To find a suitable expression for the ground-state wavefunction tp g we use the results of the preceding 
section, i.e., consider the case where the BEC is confined by a periodic potential of the form 

V(x) = ~V a cos(2k L x) (17) 

(we choose the factor of 2k^ since the potential created by an optical lattice of wavevector fc^ would 
create such a potential JhJ). In addition, we consider a very weak probe beam and neglect the four- 
wave mixing effect due to the interference between the probe laser and the lattice lasers 111]]. As a 
result, the ground state of the BEC can effectively be described by the Gross-Pitaevskii equation 

^ 9= {^ + V ^}^+9 gg M^g , (18) 

where p is the chemical potential. In the experimentally realized dilute Bose condensates the inter- 
action energy of the two-body collisions between atoms is usually large, so that one can perform the 
Thomas-Fermi approximation by neglecting the kinetic energy. This transforms Eq. (|l^) into a simple 
algebraic equation whose solution is of the form |?/j g (:E)| 2 = po + Pi cos(2ki J x). This solution is valid 
for all x if the optical potential is not too strong, so that p\ is smaller than po- If Pi is n ot too close 
to po we furthermore can simplify the wavefunction to 

^ g {x) ^ yfp~o~+ -^—cos{2k L x) (19) 
2 VPo 



(this time-independent expression is valid in a frame rotating at frequency E g /%). We remark that 
this expression produces qualitatively correct results for the lowest photonic band even if the kinetic 
energy is not negligible or p\ k, p Q because the corresponding correction essentially introduce higher 
coefficients in the Fourier series of Eq. (|l9|). Since these higher coefficients do only couple higher 
bands, they do not affect the results for the lowest band. 

Since ip g does not depend on y and z it is advantageous to rescale the wavefunctions as if) — * Lj_tjj, 
where L± is the typical extension of the BEC in the y and z direction. This guarantees that the 
one-dimensional integral J dx\ip\ 2 is dimensionless and can be interpreted as a particle number. In the 
actual calculations this rescaling leads to the appearance of various factors of L±. L± will not enter 
the final results, though. 

Introducing the (classical) field f^ p )(a;) := d-ujA^p\x) for the probe laser's Rabi frequency Eqs. (Q) 



and (|9j) can be reduced to the polariton equations of motions 



ie_pr 


ihlPe = + ^rcs + Jit l^l'j V>e ~ ity (20) 

lhtt P l = fl(j^ + l^^Ve} , (21) 

where tpg (divided by L±) is given by Eq. (^). Since the density of excited atoms should be very 
small we have neglected two-body collisions between excited atoms (g ee = in Eq. (]?])). Because tf> g 
only depends on x, and since we consider the case that ipe and do not depend on y and z either, 
the transverse delta function of Eq. (^) can be reduced to an ordinary delta function (to prove this 
one can transform Eq. (|9|) to momentum space). 

For later use it will be convenient to consider the solutions {x\4>) :— (ip e (x), Cl^(x) ) of Eqs. ( |2(i| ) 
and @ as elements of a polariton Hilbert space with conserved scalar product 

dx L' e *^ e + ^ll n (P)'^-^(P)\ (22 ) 
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Physically the quantity (<fi\(j>) is related to the number of excitations (number of excited atoms plus 
number of photons) in our system. It is a conserved quantity because of the rotating-wave approxi- 
mation made in Sec. II. Eqs. ( p0[ ) and ( f2l| ) can be rewritten in the form ihdt\4>) = H po \\<j)) with the 
polariton Hamiltonian 

tfpoi := \(1 + ^) + T^res + |f l^l 2 } + ^(1 - ^)hu - iH g a + + ^L^* g a- , (23) 

where &i are the Pauli matrices. We remark that H po \ is Hermitcan with respect to the scalar product 



©, i.e., {ct>'\H poi <t>) = (ffpoi^l^. 

To derive the polariton band-structure we have to find the eigenvalues of the operator -ff po i. Since 
tp g is periodic -ff po i commutes with the operator of discrete translations of amount 7r/fci and thus has a 
common set of eigenvectors with this operator. The eigenvectors \4> n ,q) therefore can be characterized 
by two quantum numbers n £ {Q, 1, 2, •••} and q £ [— kL,ki] which denote the band index and 
the quasi-momentum, respectively. The eigenvalues of the discrete translation operator are given 
by exp[ig7r/A;L] and belong to eigenvectors which are simply given by momentum eigenstates with 
momentum Hk m := h{q + 2mki) for integer m. The eigenvalues Tiu) n ^ q of the Hamiltonian can be 



found by expanding Eqs. ( |20[ ) and (|2l]) in this basis and searching for stationary solutions. The 
corresponding equations, 



Hu;ip e (k m ) = I + fibres j V'e(fcm) + j§- |poV>e(fcro) + y[^ e (fe m+1 ) + VeC^m-l)]} 



-inL ± ^rt p Xk m ) - inL ± -£=[rt p \k m+1 ) + r»( p )(fc m _!)] (24) 



hurt p \k m ) = hc\k m \rt p \k m ) + i ld ^ 2c } kml fv^V-e(*m) + ^=lMk m+ i)+Mk m -i)] \ , (25) 



can easily be solved numerically. Fig. lc) shows the resulting band structure near the upper band 
edge of the lowest frequency band for a condensate of density po — l.lpi = 10 14 cm" 3 . In order to 
describe the limit of a photonic instead of a polariton band-structure we have assumed a very large 
detuning of Al — 100 GHz of the lattice. However, the results given below do not change very much 
if a smaller detuning is assumed. We furthermore have set |dj_| ~ eao, with e being the electron's 
charge and ao denoting Bohr's radius. The wavevector of the lattice was taken to be k^ = 10 7 m" 1 . 

An excellent analytical approximation for the band structure can be made by assuming that for 
q £ (0, k L ) only the modes ^ p \q), ^ p \q — 2k L ), tp e (q), and ifj e {q~2ki) are important. The problem 
is then reduced to finding the eigenvalues of a 4 x 4 matrix. For q — k^ the eigenvalues have a simple 
form and allow to derive the following expression for the band gap Aw separating the two lowest 
energy bands, 

Auj = s + - s- , (26) 

where we have defined the frequencies vt :— \d±\ 2 pi/(2heo) and introduced the abbreviations s± := 
^/(Al/2) 2 + id ICS (y/Uo ± y/Ui/A) 2 . For a large detuning |A_l|, i.e., in the limit of a photonic band gap, 
this expression simplifies to Alj — ujoi>i/\Al\. For the numerical values given above the band gap 
takes the value Aui w 40 GHz. 

For q ^ k^ the band structure is given by a rather complicated expression. We therefore have 
further simplified the analytical result by fitting it to a square root ]l5| ]. The lowest polariton band 
then takes the form 



^0,q ~ Wo,max + V- y/c 2 (\q\ - k L ) 2 + V 2 , (27) 

where 

|Ai| . . 

Wo,max — w rcs ^ s + l^oj 

denotes the upper edge of the lowest frequency band, and 

Auj 



s+(s + + s-) (29) 



V^O, max — w resJ 

determines the curvature of the band. In the limit of a large detuning this simplifies to v ~ Aw/2. 
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V. THEORY OF LOCALIZED DEFECTS 



In the previous section we have studied polariton band gaps of light generated by the lattice con- 
densate in its ground state. However, in general the BEC might be in a state corresponding to a 
(coherent) elementary excitation which usually are not periodic. Thus, we expect defects in the lat- 
tice condensate. As is well-known from solid state theory a defect or an impurity in an otherwise 
periodic potential can lead to defect states, i.e., states whose energy eigenvalue lies inside the gap be- 
tween two energy bands (see, e.g., Ref. |Ls|] ). In the system under consideration this phenomenon could 
be exploited to acquire knowledge about non-periodic elementary excitations of the Bose condensate 
by observing light propagation through the BEC. In addition, defect theory can also be applied to 
study the back-reaction of an excited BEC on the optical potential confining it. 

The existence of defect states for photonic band gaps has been examined in the microwave regime 
for ordinary dielectric materials [ p^lp^ ] . The method of calculation that we adopt is closely related to 
the Green's function approach of Ref. ]l6| . 

Specifically we consider the situation that the condensate's wavefunction is given by ipa(x) — 



ipg(x)+5ip g (x), where ip g (divided by L±) is given by Eq. (19) and 8ip g describes a coherent elementary 
excitation of the condensate which we assume to be localized in the x-direction. This allows us to 
estimate the resulting energy eigenstates by adapting the Koster-Slater model to the case of 
polariton band gaps. 

Our aim is to find solutions of the equation 

hu de { cct \(t)} = (-ffpol + #defcct)|0) , (30) 

where 

^defect := -ihSip g a + + ^ rfj l 2 u!§ip*a- (31) 

describes the influence of the elementary excitation Sip g and ^defect is the defect eigenfrequency. 
To find the solutions of Eq. (|3^) we expand \<p) in terms of Wannier functions, 



/2k 



L J-k L 



dqe -™i/k L j^) t (32) 



where v is an integer number. The functions W n ^{x) are localized around the lattice point itv/kh 
which makes them a convenient tool to study localized defects. Inserting \<j>) = ^ n v <$> n ,v\Wn,v) into 
Eq. (BCj) one easily deduces the equation 



+ ^{W n , u \H Ae i ect /n\W m , li )<j) m , ll , (33) 

with 



1 



„ .- — , dqe** v -M k *w n , q . (34) 

£Kl J-k L 

Since the Wannier basis is countable Eq. (^) can be interpreted as a matrix equation. In particular, 
one now can exploit the fact that both the Wannier functions and i/dofoct are localized. This implies 
that, at least approximately, only a finite number, say an N x N submatrix, of the matrix elements 
(W n>v \H deiect \W mifl ) is nonvanishing. 
Introducing the Green's function G := (wdefectl — Hpoi/fr)^ 1 whose matrix elements are given by 

(W n ,u\G\Wm,iJ,) = S n ,m-pr7— / dq— =: Sn tm G n ,u-n (35) 

-SKL J-k L ^defect — ^n,q 

the eigenvalue problem can be reduced to 

<\>n.y = ^ GnM-n(Wn ifl \Hdcfcct/fl\Wm,,\)(j) m ,\ ■ (36) 
m,/i, A 

As is well-known in solid state theory [|l6| it is sufficient to consider only the eigenvalue problem of 
the N x N subspace where the matrix elements of -ffdefect are non-zero to derive the frequencies of 
the defect states. 
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VI. KOSTER-SLATER MODEL FOR PHOTONIC BAND GAPS 



Having derived the matrix eigenvalue equation ( pq ) for the determination of the defect frequency 
it is straightforward to apply the Koster-Slater model to the problem at hand. In this model the 
assumption is made that both the Wannier functions and the perturbation U np are localized in such 
a way that only one matrix element of the perturbation is nonzero, 



(Wn,v\Hdefect\W m ,ij.) = Uo6 n ,o8 m , 0^,0^,0 



(37) 



This model is not valid if Sip g is too strongly localized (i.e., on a scale much smaller than the lattice 
spacing 7r/fc[,) p9|| , but should produce qualitative estimates of defect frequencies for moderately 
localized perturbations. 



Inserting Eq. (37) into the eigenvalue equation (p6|) and using Eq. (R5J) to evaluate the only relevant 



matrix element of G, (Wo,o|C?| Wo,o), we find that the defect frequency Wdefect has to fulfill the condition 

i= u ^r * . (38) 

a ZKl J_k L ^defect — w 0,q 

The integral can be calculated exactly for the photonic band of Eq. (p^), but since the resulting 
expression is somewhat complicated it is more instructive to use the following approximation which 
is valid if the defect frequency ^defect is close to the upper edge lu max of the lowest frequency band, 



1 f kL dq 




2&L J-k L ^defect — ^0,q y 2(cL>dcfcct ~ ^0,max) 

Inserting this into Eq. (|38|) we find for the frequency of the defect state the expression 



(39) 



7T V Uq 

^defect — Wo,max ~ ■ (40) 



It is of interest to know how large this frequency difference is for realistic systems. To achieve this 
we first have to estimate the value of U = (W^ol-ffdefectlW^o) = ^ I dx{WnW e 8ijj* - W*Wo,8ip g } 
with (x|Wo,o) = (W e (x), Wq (x)). A rough estimate for this integral can be made by setting both 
Stjjg and the Wannier function to be constant over one wavelength and to be zero outside this 
range. The normalization condition (Wo.o|Wo.o) = 1 then leads approximately to 1 = A£{|W e | 2 + 
\Wn\ 2 / (uil\cI±\ 2 /2heoL 2 L )} . Using this condition Uq takes its maximal value for W e = X/{y2L±), 



\dj> 2 
' 2he L 2 } 



l^MVW^^TT- ( 41 ) 



Assuming a defect amplitude of Sip g — e^fp^L^ over the extent of Wo,o(x), where e is small compared 
to one, we can derive an estimate for the defect frequency (|4^) of 

^defect — Wo roax « 27T 2 £ 2 Vq . (42) 

Using the same numbers as in Sec. IV (that is, p L = 10 14 cm -3 and = 10 11 s _1 ) as well as e = 0.3 
this frequency difference can be shown to be of the order of 100 Hz. Though this number is too small to 
be measurable it should be pointed out that it applies only to a defect corresponding to an elementary 
excitation over one wavelength. A different type of defect can produce a much different result. For 
example, if we do not consider a weak elementary excitation but a strong localized excitation we can 
estimate its effect but assuming a larger value for e. A threefold increase of the local density (e = 3) 
would lead to a defect frequency ^defect — wo.mm of about 10 KHz, for instance. 
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VII. CONCLUSIONS 



In this paper we have analysed the interaction of a lattice (or "crystalized" ) Bose-Einstein conden- 
sate with largely detuned laser beams. We have derived a periodic solution of the coupled equations 
of motion, corresponding to a free BEC and a standing lattice laser beam. We found that, if the 
condensate is in its ground state, these equations decouple and the effect of the lattice laser beam on 
the condensate is not affected by the condensate itself (no back-reaction). In this situation it is thus 
equivalent to consider a condensate in some external periodic potential with the same periodicity. 

Building on this result we then assumed that the condensate is confined by an external periodic 
potential. Since the condensate's ground state is then periodic, too, it forms a kind of periodic 
dielectric. A probe laser beam propagating through this dielectric will then experience the formation 
of photonic band gaps. We have analysed this situation using the concept of polaritons, i.e., entangled 
superpositions of excited atoms and photons. 

If the condensate is not in its ground state but in a state corresponding to a localized elementary 
excitation the periodicity of the system is perturbed. This leads to the formation of defect states 
inside a polariton band gap. 
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FIG. 1. Band structure for polaritons near the band edge k^. Displayed is 10 4 (w — w rcs )/aj r0 s versus 
10 4 (g — fci)/fci. Fig. a) shows the free dispersion relation for the relevant field components in absence of 
any interaction. For a homogeneous BEC, ip g = ■s/po, a mixing of excited atoms and photons leads to the 
formation of two separate avoided crossings (solid and dashed lines in b)). The labels indicate the asymptotic 
form of the polariton modes, where they correspond to either excited atoms (tpe) or photons (f^ p ^). For a 
periodic ground-state BEC the two avoided crossings are combined to form band gaps (Fig. c) . 
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